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Abstract 

We study Einstein's equations with matter in hydrostatic equilib- 
rium in the nonstandard gauge which was recently investigated in the 
vacuum case. We obtain spherically symmetric solutions for any given 
rotation curve. These solutions can be used to describe the bulge of a 
galaxy without assuming dark matter. 
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1 Introduction 



What is nonstandard general relativity ? In standard GR one studies solu- 
tions of Einstein's equations up to equivalence classes under diffeomorphisms 
[1]. In [2] we have investigated the equivalence class of the Schwarzschild 
solution. Surprisingly enough we have found that different solutions in this 
class are physically inequivalent, that means observable quantities as the cir- 
cular velocities are different. Consequently, for physics the classification by 
means of diffeomorphisms is too coarse; a more detailed analysis is required 
and this is the matter of nonstandard GR. 

In [2] we have introduced a physical gauge fixing by requiring that an 
important observable, namely the circular velocity V(r) has prescribed val- 
ues. This leads to solutions completely different from Newtonian theory. We 
have studied the simplest situation of a static, spherically symmetric metric 

ds 2 = g^dxPdx" = e a dt 2 - e b dr 2 - r 2 e c {d$ 2 + sin 2 Mcf> 2 ) (1.1) 

in spherical coordinates (t, r, <j)) which are defined by astronomical mea- 
suring techniques and we shall keep these coordinates throughout. The 
functions a, 6, c depend on r only, and we take c(r) / which is the es- 
sential difference to standard general relativity [3]. We have studied the 
vacuum solutions in [2] . Using them we have constructed disk- like solutions 
with a somewhat unphysical energy-momentum tensor. Here we consider 
the simplest physical energy-momentum tensor 

# = diag(e,-p,-p,-p), (1.2) 

which describes ordinary matter in hydrostatic equilibrium. Since ds 2 in 
(1.1) is assumed diagonal it makes no sense to consider more complicated 
t^ u here. 

In the next section we set up the corresponding field equations, and, 
taking energy-momentum conservation into account we find the same kind 
of algebraic dependence of the equations as in the vacuum case [2]: one of 
the three functions a, b, c remains undetermined. We live with this and take 
again the circular velocity V{r) as being given, to determine the metric com- 
pletely. There results two coupled first order ordinary differential equations 
for a metric function and the pressure p(r). In standard general relativity 
we have only one equation for p(r) instead ([1], Sect. 11/1). In Sect. 3 we 
solve these equations for a spherical bulge of a galaxy where the circular 
velocity is proportional to r for small r. We find the pressure which is finite 
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and continuous at r = 0. However, p{r) cannot be expanded in powers of r 
at r = because it involves the integral logarithm. In the last section we 
consider a slightly more general energy-momentum tensor with dust. This 
gas plus dust problem has no solution in the static case. 



2 Einstein's equations in a spherically symmetric 
nonstandard gauge 

The Einstein tensor for the metric (1.1) has been calculated in [2]. Inserting 
the energy-momentum tensor on the right-hand side we obtain the following 
three field equations 



Gtt = e 



a—b 



-c" - -c' 2 + h'c' + -(&' - 3c')} +\(e a - c - e a - b ) = «t„ (2.1) 



r J V- 
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G rr = \a V + V + cO + ^ + ^(l-e b - c ) = 



G m = ye c - fe [a" + c" ~(b' -a' -2c') + ^a' 2 -a'b' + a'c' -b'c' + c' 2 ) 



1 



(2.2) 
(2.3) 

where k = 87rG/c 4 and G is Newton's constant. The prime always means 
d/dr. By combining (2.1) with (2.3) we write these equations in the form 

c" = -\d 2 + h'c' + -(&' - 3c')} +\{e b ~ c - 1) - Kt tt e b - a (2.4) 



c' 2 a 



11 Ik 

a" = ^-^ a '-b'+c')+-(c'-a')+-(l-e b - c ) + ^t m e b - c +Kt tt c b - a . (2.5) 

As in [2] we eliminate b(r) using (2.2) 



/ 1 a'c' c' 2 a! + c' \ 

6 = c + 21ogr + log — + — + — + Kt rr ) (2.6) 

\r z 2 4 r ' 



and 



e b ~ c = l + r 2 + +r{a' + c') - nr 2 t r 



M 



(2.7) 



b ' = C ' + r + lM 



a"(2r 2 c' + 4r) + c"2r z (a' + c' + -)- 4(a' + c) - - - A K r z t 



2.1 



(2.8) 
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Then from (2.5) we get 



a' + d + 2/r 



( , 2 4c' c' 2 4 
c' + - + — + — + 



K . 

, , ; v 4 — t rr ) + 

r or a ffr a 



2M 2 
+^- + ^t + 2k 



i' 



4k 



M 



r 2 a' + c' + 2/r a'r 2 a' + d + 2/r 
Substituting this into (2.8) we find 



(2.9) 



a 4 k r / 2 \ 

b 1 = 2— + a' + 2d + - - 2— i rr + M (-^t^ + t tt e c_a 
a' r a' L \r z ' 

This allows to simplify (2.9) as follows 



„ a" , , 2 X d 2 a' + d , , 4 

c" = — (c + - + — + 2 + a'c + ^ + nt 5 

a' r 2 r r z 



where 
and 



-t r r~\ 



2tL 



AM 



a' + d + 2/ra' a' + d + 2/r a'r 2 (a' + d + 2/r) 



(2.10) 



(2.11) 



U (2.12) 



U = t rr + M^t m + tae C ' a ^j . 

Next we eliminate b in (2.4) and find a result similar to (2.11): 



2, d 2 „a' + d , , 4 



+ a c H — ^ — «i c 



(2.13) 



(2.14) 



with 



12 12 
t c = t tt e c - a M [l + - (d + -)] +t rr [l + - (d + -) 



+ 



2M / , 2x 

+W-r-r( c + -)■ 

a r z V r/ 



(2.15) 



To see whether this is identical with (2.11-13) we insert the components of 
the energy-momentum tensor (1.2). Lowering one index in (1.2) we have 

t tt = e a g, t rr = e b p, t m = r 2 e c p. (2.16) 

Then t c (2.15) is equal to 

fe = e 'M 1 +^'+7)M 1 + 5( c ' + ;)]}- (2 - 17) 
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On the other hand t& (2.12) contains the derivative of t rr . This is the point 
where energy-momentum conservation comes in. In hydrostatic equilibrium 
we have ([1], eq.(5.4.5)) 



1/2 

d r p = ~{p+ g)d r log g t l , 



which gives 



p'(r) = -(p + q) — . 



(2.18) 
(2.19) 



Now we are able to calculate 



t' rr = e b (p' + pb') 



= e 



~(P + s)tt +p[2— + a+ 2c + 



i-U) 
a' ) 



, . . - . , . (2.20) 

2 c \ a' r a' Vi y ' 

where (2.10) has been used. Substituting this into (2.12) we finally end up 
with the simple result 



f 

-p-g- -^(g + 3p) 



(2.21) 



Here we have introduced the function 



c' + 



(2.22) 



or 



f(r) = c(r) + 2 log ■ 



(2.23) 



which was already used in [2]; r c is a constant of integration. Comparing 
(2.21) with (2.17) we see that t$ = —t c . Consequently, the two field equations 
(2.11) and (2.14) are identical which is the same degeneracy as in the vacuum 
case [2]. 

As in the vacuum case we lift the degeneracy by requiring that the cir- 
cular velocity V(r) has given values. The latter is equal to ([2], eq. (2.12)) 



V 2 (r) = u = 



(2.23) 



This enables us to eliminate a{r) in favor of u(r) and f(r). From (2.19) we 
obtain the differential equation for the pressure 



u 



P 



(P + Q)f 



(2.24) 
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Next we rewrite (2.14) in terms of / and u. It is a remarkable property 
that the second derivatives drop out, so that we get a first order differential 
equation for /: 

f,2 <- ' ^ ' *' U ' -nt c . (2.25) 



This equation is coupled to (2.24) because t c contains the pressure p. In 
standard general relativity on has only one first order differential equation 
for the pressure ([1], eq.(11.1.13)). 

To work out t c in detail we start from 



e b = rV 



( r 2 



a'c> c' 2 



a' + c' 



Inserting t rr = p exp b and 



we can solve for 



This finally leads to 



e b = 



■f' 2 (- + -) 

l + Kr 2 efp J V2 4/ 



r 2 c ef 



(2.26) 



P 



( 1 + -) 

V II./ 



3\ (, 1\ 



l 



4^1 + Kr 2 e->p 



(2.27) 



Substituting this into (2.25) gives us a closed system of first order differential 
equations for f(r) and p(r). As in standard general relativity an equation 
of state q = g(p) must also be known. 



3 Integration of the equations for small r 

Regarding the application of our spherically symmetric metric to real galax- 
ies it is a pity that spherical galaxies are rare. If one analyzes disk galaxies, 
for example M33 [4] [5], one represents the stars in the disk by some den- 
sity distribution which represents the measured surface brightness. However, 
then quite often there remains an emission excess in the central region. This 
is attributed to a nucleus or bulge, and it is not bad to describe this bulge 
by a nonstandard spherically symmetric metric and corresponding matter 
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density and pressure. Therefore we want to solve our above equations for 
small r. 

The behavior of the rotation velocity V{r) for small r is well known [6], 
V(r) vanishes linearly for r — > 0, therefore we assume 

u{r) =u 2 r 2 + 0(r 3 ) (3.1) 

where u 2 is a positive constant. Regarding the matter density g we set 
g = const throughout to find the most simple solution. Then as the zeroth 
approximation we also put the pressure p to be constant p = po, but later on 
we calculate corrections to this. Under these assumptions we find a singular 
behavior for t c (2.27) 

tc = % (3-2) 

where t 2 is another positive constant. Now equation (2.25) assumes the 
simple form 

r 2 /' 2 +4rf = 2Kt 2 e f . (3.3) 
We solve the quadratic equation 



r-f = -2±^4 + 2Kt 2 ef (3.4) 



and we integrate this equation by separating the variables 

dr_ df 

r ~ -2(1 T Jl + %t 2 ef' 



(3.5) 



Only the upper minus sign in (3.5) leads to a physical solution. After 
integration we arrive at the equation 



r 1 



1 + 1/ 



1 



-21og- = -log-^+ (3.6) 

77 2 1-y y-1 



where 



1 + Jt 2 e/ « 1 + -Met (3.7) 
2 4 

and Tf is a constant of integration. We shall soon see that exp/ is small 
for small r such that we can expand the square root in (3.7). For the same 
reason we can neglect the logarithmic term in (3.6) so that 

_21 og - = -L e -/ (3. 8 ) 

Vf Kt 2 



7 



or 



J _ 



(3.9) 
nt2 log \r/rf\ 

for r < rf. Indeed this goes to zero for r — > 0. Similarly from (3.4) we find 



f = -(-l + Jl + -t 2 ef)*-t 2 - 



r log \r/rf\ 



We also note 



f(r) = log — - log log 

Kl2 



which follows from (3.9). 

Next we must integrate the equation (2.24) for the pressure 



Using 



r 2 f' = 



log \r/rj 

the solution of the homogeneous equation is given by 



p(r) = pi exp 



4 J ^r f 

where p\ is a constant of integration and 

z 

dx 



\iz 



I 



log a; 



(3.10) 
(3.11) 



(3.12) 



(3.13) 



is the integral logarithm. Then the solution p(r) of (3.12) with p(0) = po is 
equal to 



V = (Po + 0) exp 



(3.14) 



The integral logarithm has no expansion in powers of r at r = 0, so we leave 
(3.14) as its stands. Still the pressure is continuous and differentiable at 
r = 0. 
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We also want to calculate the metric, although it does not have a direct 
physical meaning. From (3.9) we immediately find 



e c = m 2 ef = 2Tc (3.15) 

V r / Kt 2 r 2 \ log r/rj\ 

The it-component follows from 

a =uf = n — r~\ 
| log r/rf | 

which gives 

e a = ^ a exp(^h(-) 2 ). (3.16) 

Finally, exp6 is obtained from (2.26). In leading order we have 

r 2 1 

e b = -^—5- (3.17) 

2k^2 r z \ \ogr/rf\ 6 

So the metric has a singularity at r = 0, but there is no horizon at some 
finite r as in the Schwarzschild solution. 



4 The gas plus dust problem 

The energy-momentum tensor (1.2) which we have considered so far can be 
interpreted as coming from a gas in hydrostatic equilibrium. Now we want 
to add dust with density £o to this gas so that we assume a slightly more 
general 

t£ = diag(e + 0o,-p,-p,-p). (4.1) 

By definition the dust has no pressure; that means p is the gas pressure as 
before. The additional £o violates energy-momentum conservation, therefore 
we expect some obstruction. 

It is not hard to repeat the calculations in sect. 2 with this more general 
energy-momentum tensor (4.1). Instead of (2.17) we have 

t c = e b {( 6 + Qo)(l + £)+p(l + 3^)}, (4.2) 

but the energy-momentum conservation (2.19) remains unchanged because 
it holds for the gas only. Regarding t$, equation (2.12) still holds, but (2.13) 
becomes 

*4 = e b (g + Q + 3p). (4.3) 
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which causes a change in t' rr (2.20). Substituting all this into (2.12) we 
finally obtain 

*5 = -tc + Qve b ^— r (4.4) 

Here the last term destroys the previous degeneracy in the Einstein's equa- 
tions. Moreover one easily sees that there is no solution of the equations for 

eo^o. 

This negative result has a clear physical interpretation: static dust is 
impossible. Dust is also used to represent stars in a galaxy model. Then 
the physics is pretty clear: the stars must move. This gives additional 
components in the energy-momentum tensor (4.1) and, as a consequence, 
the spherical symmetry is destroyed. 
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